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Abstract 


In  this  report  we  show  how  magnetic  noise  occurring  on  an  ocean  floor  magnetometer  can 
be  reduced  using  a  collocated  pressure  sensor  and  a  separate  more  distant  magnetometer. 
To  achieve  this,  a  multi-channel  coherence  technique  is  applied  to  the  measurement  data. 
The  removal  of  magnetic  noise  on  ocean  floor  magnetometers  has  several  applications,  in¬ 
cluding  mine  warfare  systems  and  ocean  surveillance  systems.  Decreasing  the  background 
magnetic  noise  by  10-20  dB  will  allow  the  detection  of  a  vessel  magnetic  field  with  am¬ 
plitude  the  same  order  of  magnitude  as  the  background  field,  and  with  frequency  content 
similar  to  that  of  the  background  field. 


Resume 


Dans  le  present  rapport,  nous  montrons  de  quelle  fagon  le  bruit  magnetique  enregistre  a 
un  magnetometre  pose  sur  le  plancher  oceanique  peut  etre  reduit  au  moyen  d’un  capteur 
de  pression  co-implante  et  d’un  autre  magnetometre,  situe  a  une  certaine  distance.  A  cette 
fin,  on  applique  une  technique  de  coherence  multicanal  aux  donnees  de  mesure.  La  sup¬ 
pression  du  bruit  magnetique  a  des  magnetometres  sur  le  plancher  oceanique  a  plusieurs 
applications,  notamment  en  ce  qui  concerne  les  systemes  de  guerre  des  mines  et  de  sur¬ 
veillance  des  oceans.  La  diminution  de  10  a  20  dB  du  bruit  magnetique  de  fond  permettra 
la  detection  du  champ  magnetique  d’un  navire  ayant  une  amplitude  du  meme  ordre  de  gran¬ 
deur  que  le  bruit  de  fond  et  un  contenu  en  frequence  similaire  a  celui  du  champ  de  fond. 
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Executive  summary 


Noise  Removal  Using  Multi-Channel  Coherence 

C.E.  Lucas,  R.  Otnes;  DRDC  Atlantic  TM  201 0-302;  Defence  R&D  Canada  - 
Atlantic;  December  201 0. 

Background:  The  removal  of  background  noise  from  bottom-mounted  magnetic  sensors 
has  applications  in  both  mine  warfare  and  target  surveillance  systems.  For  a  particular  mag¬ 
netic  sensor,  we  can  reduce  geomagnetic  and  swell  induced  magnetic  noise  by  using  data 
from  both  a  collocated  pressure  sensor,  and  a  separate  magnetic  sensor  situated  approx¬ 
imately  100  m  away.  We  use  a  multi-channel  coherence  technique  to  perform  the  noise 
removal. 

Principal  results:  For  frequencies  below  200  mHz  we  show  that  by  using  the  multi¬ 
channel  coherence  technique  we  can  reduce  background  magnetic  noise  by  up  to  20  dB. 
Using  a  remote  magnetometer  as  the  first  of  two  reference  signals,  we  reduced  the  coherent 
background  geomagnetic  noise  in  the  5  mHz-60  mHz  frequency  band.  As  a  second  refer¬ 
ence  signal  we  used  a  collocated  pressure  sensor  and  reduced  the  swell-induced  magnetic 
noise  in  the  60  mHz-200  mHz  frequency  band. 

Significance  of  results:  Decreasing  the  background  noise  on  a  magnetic  surveillance  sen¬ 
sor  by  10-20  dB  in  the  5  mHz-200  mHz  frequency  band  will  allow  target  detection  algo¬ 
rithms  to  detect  surface  and  sub-surface  targets  whose  magnetic  signatures  are  typically  of 
the  same  amplitude  and  frequency  as  the  background  noise. 

Future  work:  A  real-time  algorithm  for  the  implementation  of  the  multi-channel  coher¬ 
ence  technique  should  be  investigated.  This  technique  should  be  tested  under  different 
ocean  surface  conditions  and  geomagnetic  noise  conditions.  Determining  the  optimal  sen¬ 
sor  separations  that  maximize  the  noise  reduction  using  this  technique  also  needs  to  be 
studied. 
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Noise  Removal  Using  Multi-Channel  Coherence 

C.E.  Lucas,  R.  Otnes ;  DRDC  Atlantic  TM  201 0-302 ;  R  &  D  pour  la  defense 
Canada  -  Atlantique ;  decembre  2010. 

Introduction  :  La  suppression  du  bruit  de  fond  a  des  capteurs  magnetiques  poses  sur  le 
fond  a  des  applications  dans  les  systemes  de  guerre  des  mines  et  de  surveillance  des  cibles. 
Dans  le  cas  d’un  capteur  magnetique  particulier,  nous  pouvons  reduire  le  bruit  magnetique 
induit  par  la  houle  et  geomagnetique  a  l’aide  de  donnees  d’un  capteur  de  pression  co- 
implante  et  d’un  capteur  magnetique  distinct  situe  a  environ  100  m.  Nous  utilisons  la  tech¬ 
nique  de  la  coherence  multicanal  pour  effectuer  la  suppression  du  bruit. 

Resultats  :  Dans  le  cas  des  frequences  au-dessous  de  200  mHz,  nous  montrons  la  tech¬ 
nique  de  la  coherence  multicanal  nous  permet  de  reduire  le  bruit  magnetique  de  fond  de 
jusqu’a  20  dB.  En  nous  servant  du  signal  d’un  magnetometre  eloigne  comme  premier  de 
deux  signaux  de  reference,  nous  avons  reduit  le  bruit  geomagnetique  de  fond  coherent  dans 
la  bande  de  frequences  5  mHz-60  mHz.  Comme  second  signal  de  reference,  nous  avons  uti¬ 
lise  un  capteur  de  pression  co-implante  et  reduit  le  bruit  magnetique  induit  par  la  houle  dans 
la  bande  de  frequences  60  mHz-200  mHz. 

Portee  :  La  diminution  du  bruit  de  fond  de  10  a  20  dB  dans  la  bande  de  frequences  5  mHz- 
200  mHz  a  un  capteur  de  surveillance  magnetique  rendra  possible  la  detection,  a  l’aide 
d’algorithmes  de  detection  de  cibles,  de  cibles  de  surface  et  sous-marines  dont  la  signature 
magnetique  est  normalement  de  la  meme  amplitude  et  de  la  meme  frequence  que  le  bruit 
de  fond. 

Recherches  futures  :  II  faudrait  explorer  un  algorithme  en  temps  reel  en  vue  de  la  mise 
en  oeuvre  de  la  technique  de  la  coherence  multicanal.  Cette  technique  devrait  etre  mise  a 
l’essai  dans  differentes  conditions  de  la  surface  de  1’ ocean  et  du  bruit  geomagnetique.  II 
faut  aussi  etudier  les  moyens  de  determiner  la  distance  optimale  de  separation  des  capteurs 
pour  optimiser  la  reduction  du  bruit  a  l’aide  de  cette  technique. 


IV 


DRDC  Atlantic  TM  2010-302 


Table  of  contents 

Abstract .  i 

Resume .  i 

Executive  summary .  iii 

Sommaire .  iv 

Table  of  contents  .  v 

List  of  figures .  vi 

1  Introduction .  1 

2  Background .  2 

3  Multi-Channel  Coherence .  3 

3.1  Mathematical  Development .  4 

4  Analysis  of  Experimental  Data .  6 

5  Conclusions .  13 

References .  15 


DRDC  Atlantic  TM  201 0-302  v 


List  of  figures 


Figure  1 :  Signal  flow  diagram  for  the  mathematical  model  relating  the  test  signal 

Y  to  two  input  signals  X\  and  X2 .  3 

Figure  2:  Norwegian  Herdla  range  sensor  positions .  5 

Figure  3:  Coherence  between  BzE(t)  and  the  pressure  signal  PE(t),  East  sensor 

group .  7 

Figure  4:  Coherence  between  vertical  magnetic  fields  BzE(t )  and  Bzm^ ) .  8 

Figure  5:  Multi-channel  coherence  between  BzE(t)  and  both  7J/-;(0 .  9 

Figure  6:  Original  time  series  y(t)  —  BzE(t)  and  the  calculated  residual  signal  e(l).  10 

Figure  7:  Power  spectral  densities  Syy(f)  and  See(f)  of  the  signals  y(t)  —  BzE(t ) 

and  e(t) .  10 

Figure  8:  The  degree  of  cancellation  (dB)  between  y(t)  =  BzE(t )  and  both  signals 

Xi  (t)=BzM{t),X2{t)=PE{t) .  11 


VI 


DRDC  Atlantic  TM  2010-302 


1  Introduction 


Ocean  bottom-mounted  magnetometers  are  used  extensively  in  oceanic  surveillance  sys¬ 
tems  [1],  These  systems  are  mostly  deployed  in  water  depths  no  more  than  a  few  hundred 
meters.  Surveillance  systems  typically  consist  of  an  array  of  sensors,  or  multi-sensor  plat¬ 
forms,  spaced  according  to  the  particular  application,  such  as  port  access  or  deep  water 
choke  point  surveillance.  The  underwater  sensors  or  platforms  are  usually  cabled  to  each 
other  and/or  cabled  to  shore.  In  such  connected  arrays  the  data  from  all  sensors  is  available 
at  each  sensor  or  onshore  immediately.  Signal  processing  methods  can  take  advantage  of 
having  all  the  sensor  data  from  all  channels  available,  and  can  be  used  to  reduce  unwanted 
noise  from  particular  channels  of  interest. 

Most  surface  and  sub-surface  ocean  vessels  have  both  permanent  and  induced  magneti¬ 
zation  within  their  structures.  This  magnetization  produces  an  associated  magnetic  field, 
often  referred  to  as  a  signature,  that  can  be  sensed  by  modem  magnetometers  at  hundreds 
of  meters  range.  When  a  vessel  transits  near  a  magnetometer,  the  magnetometer  measures 
the  vessel’s  quasi-static  magnetic  field,  which  varies  only  with  time  as  the  vessel  moves 
relative  to  the  fixed  sensor. 

Underwater  surveillance  systems  that  employ  magnetometers  use  signal  processing  algo¬ 
rithms  to  detect  magnetic  signatures  from  nearby  vessels.  Reducing  the  effective  back¬ 
ground  noise  on  the  magnetometer  measurements  will  allow  the  algorithms  to  detect  smaller 
target  signatures  that  otherwise  would  have  be  hidden  in  the  background  noise. 

The  magnetic  signature  produced  from  a  transiting  vessel  measured  by  a  bottom-mounted 
magnetometer  lies  predominantly  in  the  frequency  band  from  approximately  5  mHz  to 
200  mHz.  Within  this  frequency  band  there  are  two  dominant  sources  of  background  mag¬ 
netic  noise,  geomagnetic  and  ocean  swell  induced  [2,  3].  Geomagnetic  noise  is  due  to 
electrical  currents  deep  within  the  Earth,  as  well  as  interactions  between  the  solar  wind 
and  the  geomagnetic  field.  It  is  well  known  that  geomagnetic  noise  has  high  spatial  co¬ 
herence  over  many  hundreds  of  meters  and  even  kilometers  [4] .  We  can  take  advantage  of 
the  coherence  if  we  have  multiple  magnetic  sensors  with  sensor  separations  on  the  order  of 
100  m  or  more,  with  one  sensor  in  the  same  vicinity  of  the  vessel. 

Ocean  swell  induced  magnetic  noise  is  due  to  eddy  currents  arising  from  the  movement  of 
sea  water,  and  the  associated  charged  particles,  in  the  Earth’s  magnetic  field  during  swell 
activity.  These  electrical  currents  have  a  corresponding  magnetic  field  associated  with 
them  which  contributes  to  the  noise  measured  on  ocean-bottom  magnetometers.  It  has 
been  established  that  there  is  a  physical  relationship  between  the  swell  induced  magnetic 
noise  and  the  local  hydrostatic  pressure  at  the  sensor  location  [5]  [6]. 

In  this  report  we  describe  the  use  of  a  multi-channel  coherence  technique  to  reduce  back¬ 
ground  noise  on  the  vertical  magnetic  field  measured  from  a  bottom-mounted  magnetome- 
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ter  (our  test  signal).  We  use  measurements  of  the  local  swell-induced  pressure  signal,  and 
the  magnetic  field  from  a  separate  magnetometer  approximately  100  m  away,  as  additional 
input  data  for  the  multi-channel  coherence  technique.  By  using  these  signals  we  were  able 
to  reduce  the  magnetic  noise  on  our  test  signal  by  up  to  20  dB  in  the  5  mHz-200  mHz 
frequency  band. 

2  Background 


Consider  a  vessel  transiting  within  measurement  range  of  a  bottom-mounted  magnetometer 
in  water  at  depth  on  the  order  of  100  m.  We  consider  this  sensor  to  be  our  test  sensor,  and 
the  measured  vertical  magnetic  field  our  test  signal.  The  measured  magnetic  field  is  the 
sum  of  the  target  magnetic  signature,  the  geomagnetic  background,  and  the  swell-induced 
magnetic  noise. 

To  reduce  the  background  noise  on  the  test  signal,  we  use  the  local  hydrostatic  pressure 
signal,  and  the  vertical  magnetic  field  measured  at  a  distant  magnetometer  approximately 
100  m  away.  Using  the  multi-channel  coherence  technique,  we  can  remove  from  our  test 
signal  all  signal  components  that  are  coherent  with  both  the  collocated  pressure  signal  and 
the  distant  magnetometer  signal.  We  will  refer  to  the  local  pressure  and  distant  magnetome¬ 
ter  signals  as  our  reference  signals.  Ideally  the  reference  signals  will  only  have  coherence 
with  the  unwanted  background  noise  on  the  test  signal. 

To  apply  this  technique  to  the  detection  of  a  vessel’s  magnetic  field,  we  take  advantage 
of  the  field’s  rapid  decay  with  range  R.  The  magnetic  field  of  a  vessel  can  be  accurately 
represented  as  that  produced  by  one  or  more  magnetic  dipoles.  The  resulting  magnetic  field 
amplitude  decreases  with  range  at  least  as  fast  as  l/R3.  For  a  vessel  transiting  near  our 
test  sensor,  the  associated  vessel  magnetic  field  100  m  or  more  away  will  be  significantly 
reduced,  possibly  completely  hidden  in  the  noise  or  below  the  sensor  resolution,  because 
of  the  decay  with  range.  If  this  is  the  case,  the  remote  sensor  signal  will  consist  mostly  of 
geomagnetic  and  swell-induced  magnetic  noise.  If  we  have  access  to  the  remote  magnetic 
signal,  we  can  use  the  multi-channel  coherence  technique  to  remove  from  our  test  signal 
those  components  that  are  coherent  with  the  remote  signal.  Since  geomagnetic  noise  has 
high  spatial  coherence  over  large  ranges,  we  will  be  able  to  reduce  the  geomagnetic  noise 
significantly  from  our  test  signal. 

Transiting  surface  vessels  produce  a  small  pressure  signal,  which  can  be  detected  in  shal¬ 
low  water  in  calm  low-swell  conditions  by  a  pressure  sensor.  This  pressure  signal  decreases 
rapidly  with  range,  and  is  often  undetectable  even  in  low  swell  conditions.  The  local  swell 
pressure  usually  dominates  the  pressure  sensor  signal.  Additionally,  the  pressure  signature 
of  a  vessel  typically  has  low  coherence  with  it’s  magnetic  signature.  If  we  assume  our  pres¬ 
sure  sensor  is  at  a  depth  where  the  pressure  signature  of  the  vessel  is  small  or  incoherent 
with  the  it’s  magnetic  signature,  we  can  then  safely  use  the  measured  pressure  to  remove 
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the  coherent  swell-induced  magnetic  signal  from  our  magnetic  test  signal.  As  described 
in  [5]  [6],  the  swell-induced  magnetic  noise  is  coherent  with  the  local  pressure  signal  pro¬ 
duced  by  the  swell.  We  use  measurements  of  the  swell  pressure,  and  take  advantage  of  its 
coherence  with  the  swell-induced  magnetic  noise,  as  in  [2],  to  remove  this  noise  from  our 
test  signal.  The  method  described  here  however  differs  from  [2]  in  that  we  use  both  the 
local  pressure  signal  and  a  remote  magnetometer,  and  use  a  multi-channel  coherence  func¬ 
tion  to  reduce  the  noise  on  the  test  signal.  Using  the  multi-channel  coherence  technique 
we  simultaneously  remove  from  our  test  signal  all  signals  that  are  coherent  with  the  local 
pressure  signal,  and  the  remote  (100  m)  magnetic  signal.  This  reduces  the  noise  on  the  test 
signal  and  makes  the  vessel  magnetic  signature  easier  to  detect  using  a  detection  algorithm. 

3  Multi-Channel  Coherence 


A  multi-channel  coherence  function  can  be  used  to  remove  signal  components  from  a  test 
signal  that  are  coherent  with  multiple  other  signals  [7].  Figure  1  shows  the  signal  flow 
diagram  for  the  mathematical  model  relating  a  test  channel  Y  and  two  separate  signals 
X\  and  Xi  that  are  coherent  with  Y .  In  the  signal  flow  diagram  capital  letters  denote  the 
frequency  domain  representation  of  the  signals.  The  functions  H\  and  Hi  represent  the 
linear  system  transfer  functions  that  relate  the  two  input  signals  X\ ,  X2  to  the  output  signal 
Y .  As  seen  in  Figure  1  the  signal  V  is  the  sum  of  the  outputs  of  the  two  transfer  functions, 
and  signal  Y  is  the  sum  of  signals  V  and  E.  The  residual  signal  E  is  that  part  of  signal  Y 
that  is  not  coherent  with  X\,  X2,  and  signal  V  is  that  part  of  signal  Y  that  is  coherent  with 
XUX2. 


Figure  1 :  Signal  flow  diagram  for  the  mathematical  model  relating  the  test  signal  Y  to  two 
input  signals  X\  and  X2. 

The  concept  of  application  for  this  method  is  that  signal  Y  is  the  measured  vertical  magnetic 
field  of  a  target  vessel  plus  magnetic  noise  associated  with  geomagnetism  and  ocean  swell. 
We  will  associate  signal  X\  with  the  local  pressure  signal,  and  signal  X2  with  the  vertical 
magnetic  field  from  a  reference  magnetometer.  We  assume  the  target  magnetic  signature 
at  the  reference  sensor  will  be  very  small,  or  not  coherent,  with  the  local  vessel  signature 
contained  in  Y.  The  residual  signal  E  will  contain  the  vessel  signature  contained  in  Y  with 
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all  coherent  signals  from  X\  and  X2  removed.  A  target  detection  algorithm  for  magnetic 
signatures  would  be  applied  to  the  time-domain  representation  of  signal  E. 

3.1  Mathematical  Development 

Using  cross-spectral  relationships  between  the  signal  Y  and  the  two  input  signals  Xi,  X2 
in  Figure  1,  we  can  obtain  estimates  of  the  transfer  functions  H 1,  Hi-  The  cross-spectra 
between  signals  E  and  X\ ,  E  and  X2,  and  E  and  V  are  zero  by  definition  of  the  model.  We 
denote  the  cross-spectra  as  defined  in  [7]  between  any  two  signals  X  and  Y  by  Sxy.  Taking 
the  cross-spectra  between  Y  and  the  two  input  signals  X\,  X2  we  obtain: 


Sxiy(f)  =  Sxhxi(f)Hi(f)+Sxix2(f)H2(f) 

Sx2y(f )  =  Sx2xi(f)Hi(f)+Sx2X2(f)H2(f).  (1) 

Using  the  cross-  and  auto-spectra  obtained  from  the  measured  data,  equations  1  are  solved 
for  the  unknown  transfer  functions  Hi,  H2.  Using  the  estimates  for  H\ ,  H2,  we  calculate 
the  power  spectral  density  for  signal  V  from 


Svv(f)  =  Sxix\  (/) \H\  (/)  |2  +  Sx2x2{f) \H2(f)\2  +  2Si{Sx2xi(f)Hi  (/)//£  (/)} 

=  Sx\y(f)H*(f)  +  Sx2y(f)H%(f).  (2) 


Here  *  represents  the  complex  conjugate.  Since  Syv(f)  —  Svv(  f),  and  Svv(f)  —  (Sw(f))*, 
we  can  define  the  multi-channel  coherence  function  between  Y  and  all  inputs  X,  as 

2  (f)=  \Syv(f)\2  =  (Svv(f))*Svv(f)  =  Svv(f) 

y:Xi  Syy(f)Sw(f)  Syy(f)Sw(f)  Syy(f) ' 

With  this  definition  for  the  multi-channel  coherence,  we  can  calculate  the  power  spectral 
density  of  the  residual  signal  E  from 


See{f)=Syy{f)-Svv{f)  =  (1  -y 2y:Xi(f))Syy(f).  (4) 
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Figure  2:  Norwegian  Herdla  range  sensor  positions.  Stars  denote  the  permanently  installed 
Norwegian  sensor  groups  North,  South,  East,  and  MSP.  Other  symbols  denote  sources  and 
sensors  deployed  by  DRDC  Atlantic. 


The  time  domain  residual  signal  e(t)  may  be  calculated  from  the  frequency  domain  func¬ 
tions  using  the  Inverse  Fourier  Transform  (F  '). 

e(t)=F-l{Y(f)-Xl(f)Hl(f)-X2(f)H2(f)}  (5) 

We  refer  to  the  signal  e(t)  as  the  residual  time  series  that  results  from  “cancelling”  the 
output  signal  y(t)  against  the  input  signals  x\  (t),  X2  (t). 

The  multi-channel  coherence  technique  may  be  generalized  to  N  inputs  Xi(t),i  —  1  ..N, 
allowing  further  cancellation  using  other  signals  that  have  coherence  with  the  test  signal 
y(t).  Ideally  these  other  signals  will  not  be  coherent  with  any  target  signature  contained 
in  y(t).  In  a  multi-sensor  surveillance  application  this  generalization  to  more  inputs  could 
have  a  significant  benefit  by  reducing  unwanted  noise  on  multiple  sensors  simultaneously. 
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4  Analysis  of  Experimental  Data 


Experimental  data  was  collected  at  the  Next  Generation  Autonomous  Systems  (NGAS) 
Joint  Research  Project  trial  held  in  Norway  in  September  2004.  The  trial  took  place  at  the 
Norwegian  Herdla  sensor  range  in  Hjeltefjorden,  northwest  of  Bergen. 

Figure  2  shows  the  physical  layout  of  the  Herdla  range  sensor  groups  dubbed  North,  South, 
East,  and  MSP  that  were  used  for  the  collection  of  the  data.  These  sensor  groups  are  multi¬ 
influence  platforms  that  include  tri-axial  magnetometers  and  hydrostatic  pressure  sensors. 
The  sensor  groups  are  spaced  over  100  m  apart  from  each  other,  and  their  depths  range 
from  approximately  20-85  m:  North  (49.9  m),  East  (26.7  m),  South  (84.2  m)  and  MSP 
(21.0  m).  The  data  used  in  this  analysis  was  collected  on  September  24th,  2004  and  is 
118  minutes  in  length.  During  these  measurements,  there  was  a  gale  with  wind  speeds  of 
15-20  m/s. 

For  the  data  collected  in  this  experiment  we  used  a  sample  rate  of  1 .0  Hz.  The  total  number 
of  time  samples  analyzed  per  sensor  was  NT  —  1*  1024  =  7168.  A  total  of  13  1024- 
point  data  segments  with  50%  overlap  were  used  to  estimate  the  required  spectra  using 
the  Fast  Fourier  Transform  (FFT)  and  Welch’s  Method.  Each  1024-point  data  segment 
was  windowed  using  a  Hanning  window,  then  zero-padded  to  a  length  of  NT  points  before 
taking  the  FFT.  The  resulting  spectra  were  then  averaged  over  the  13  segments  to  obtain 
the  required  auto-  and  cross-spectra.  Thus  all  spectral  estimates  in  equations  1  through  5 
including  the  transfer  functions  H\(f)  and  H2  (/)  are  NT  frequency  points  in  length.  The 
spectra  Xi(f),X2(f)  and  Y (/')  are  obtained  by  taking  the  NT -point  FFT  of  each  of  the 
corresponding  time  series. 

In  our  analysis  we  use  as  our  test  channel  the  vertical  magnetic  field  BzE(t )  measured  from 
the  Herdla  East  sensor  group.  From  this  channel  we  remove  signal  content  that  is  coherent 
with  those  signals  measured  simultaneously  from  the  East  sensor  group  pressure  sensor 
P#(f)  and  also  from  the  MSP  vertical  magnetic  sensor  BzM{t).  We  relate  our  measured 
data  with  the  signal  flow  diagram  of  Figure  1  by  assigning  y(t)  —  BzE(t),  *1  (?)  =  Pe(t), 
and  x2(t)  =  Bzm  (t )  • 
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Coherence  7*(f)  (Bz.Pres)  East  Sensor 


Figure  3:  Coherence  between  Bzt(t)  and  the  pressure  signal  PiAt  ),  East  sensor  group. 
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Coherence  y2( f)  (Bz  East,  Bz  MSP) 


Figure  4:  Coherence  between  vertical  magnetic  fields  Bze(J  )  and  Bz.m(j)- 


The  measured  data  at  the  East  sensor  group  shows  a  physical  relationship  between  the 
local  swell  induced  magnetic  noise  and  the  pressure  signal.  This  relationship  is  apparent 
in  Figure  3  where  we  have  a  high  degree  of  coherence  between  the  vertical  magnetic  field 
BzE^t)  and  the  Pressure  signal  Pe  (l )  at  the  East  sensor  group,  particularly  in  the  60  mHz  to 
200  mHz  frequency  band.  Here  we  have  used  the  coherence  function  y %x{f)  between  two 
signals  x(t),y(t),  defined  in  terms  of  their  auto-  and  cross-spectra  as 


\Sxy(f)\2 

Sxx(f)Syy(f) ' 


(6) 


In  Figure  4  we  show  the  coherence  between  the  vertical  magnetic  fields  Bze^)  and  Bz,m(j)- 
The  horizontal  distance  between  these  sensor  groups  is  approximately  100  m.  The  high 
coherence  below  60  mHz  is  due  to  the  high  spatial  coherence  of  the  geomagnetic  noise  in 
this  frequency  band  over  the  sensor  group  separation.  Geomagnetic  noise  below  60  mHz 
can  show  high  coherence  over  several  kilometers. 

In  Figure  5  we  show  the  multi-channel  coherence  between  BzEit )  and  both  Bzjvtit )  and 
Pe  ( t ),  defined  by  equation  3.  The  multi-channel  coherence  is  high  in  both  the  geomagnetic 
and  the  pressure  swell  noise  frequency  bands,  and  extends  in  frequency  up  to  200  mHz. 
Using  equation  5,  we  can  solve  for  the  residual  time  series  e(t).  Here  e(t)  is  the  East  sen¬ 
sor  group  vertical  magnetic  field  BzE(t )  with  the  coherent  noise  from  BzM(t )  and  PeO) 
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Coherence  y2( f)  Bz  East  (Bz  MSP,  P  East) 


Figure  5:  Multi-channel  coherence  between  BzEit)  and  both  Bzm (j),  PE{t). 


removed.  In  Figure  6  we  compare  the  residual  signal  e(t)  with  the  original  vertical  mag¬ 
netic  field  test  signal  y(t).  Figure  6  clearly  shows  that  a  target  signature  of  only  a  few 
nano-Tesla  can  be  easily  detected  in  the  residual  signal  e(t).  There  were  no  target  transits 
in  this  data  set. 

Figure  7  shows  the  corresponding  power  spectral  densities  Syy(f),See(f)  of  the  original 
test  channel  y(t)  —  BzE(t)  and  the  residual  signal  e(t  )  after  performing  the  cancellation 
with  the  signals  x\(t)  —  Bzwi(t )  and  *2 (t)  =  In  Figure  8  we  show  the  degree  of  can¬ 
cellation  between  y(t)  —  Bze(J)  and  both  x\(t)  =  Bz.mU)  and  X2(t)  =  PeO)-  The  degree 
of  cancellation  is  defined  to  be  the  difference  in  the  decibel  values  of  Syy(f)  and  See(f) 
of  Figure  7.  It  is  the  multiplying  factor  increase  in  SNR,  expressed  in  decibels,  result¬ 
ing  from  the  noise  removal.  Up  to  20  dB  of  noise  has  been  removed  from  the  measured 
vertical  magnetic  field  in  the  5  mHz  to  200  mHz  frequency  band.  This  frequency  band 
corresponds  to  that  associated  with  the  magnetic  signature  of  a  transiting  vessel  measured 
from  a  bottom-mounted  magnetometer.  A  target  detection  algorithm  applied  to  the  residual 
signal  e(t)  will  have  greater  detection  capability  due  to  the  factor  of  10  (20  dB)  increase 
in  SNR  .  This  applies  to  the  detection  of  both  nearby  vessels  with  small  magnetic  source 
strengths,  and  vessels  with  larger  magnetic  source  strengths  that  are  at  a  greater  distance 
away. 
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Signal  y(t)  and  Residual  e(t) 


Figure  6:  Original  time  series  y(t )  =  Bze(^)  and  the  calculated  residual  signal  e(t  )  with 
coherent  signals  from  Bzif(t),  BeU)  removed.  A  target  magnetic  signature  of  only  a  few 
nano-Tesla  could  now  be  easily  detected  on  e(t). 


Syy(f)  and  See(f)  Spectra 


Figure  7:  Power  spectral  densities  Syy(f)  and  See(f)  of  the  signals  y(t)  —  Bz.e(J)  ande{t). 
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Degree  of  Cancellation 


Figure  8:  The  degree  of  cancellation  (dB)  between  y(t)  =  Bze(j)  and  both  signals  x\  ( t )  = 
BzM(t),  X2 (t)  —  Pe  (?) .  This  measure  is  the  multiplying  factor  increase  in  SNR,  expressed 
in  decibels,  resulting  from  the  noise  removal. 
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To  optimize  the  performance  of  the  multi-channel  coherence  technique  in  reducing  noise  on 
the  test  signal,  without  also  reducing  the  target  signature  within  the  test  signal,  the  reference 
signals  must  be  chosen  carefully.  If  possible  the  reference  signals  should  be  chosen  such 
that  they  are  coherent  with  the  background  noise  on  the  test  signal,  and  not  coherent  with 
the  vessel  signature  of  interest.  This  is  likely  to  be  true  in  practical  scenarios  based  on  the 
concept  of  usage  described  here  because: 

1.  The  detection  range  of  pressure  signatures  is  known  to  be  much  smaller  than  that  of 
magnetic  signatures.  For  a  surveillance  array  deployed  at  depths  greater  than  20  m, 
the  target  pressure  signal  from  a  collocated  pressure  sensor  will  be  dominated  by 
the  local  pressure  swell.  The  pressure  signal  will  have  no  coherence  with  the  target 
magnetic  signal  on  the  test  sensor. 

2.  For  the  magnetic  signal,  our  reference  signal  would  contain  very  little  of  the  tar¬ 
get  magnetic  signature  because  of  the  rapid  fall-off  rate  with  range  of  the  magnetic 
field.  Essentially  the  target  signature  is  not  present  at  both  sensors  simultaneously. 
The  remote  reference  signal  should  then  only  consist  of  background  geomagnetic 
and  ocean-swell  induced  magnetic  noise,  and  be  incoherent  with  the  target  magnetic 
signature. 


12 


DRDC  Atlantic  TM  2010-302 


5  Conclusions 


We  have  shown  that  the  multi-channel  coherence  technique  can  be  effective  in  reducing 
unwanted  noise  on  a  test  signal  by  removing  those  signal  components  that  are  coherent 
with  other  available  signals.  For  a  test  signal  we  used  the  vertical  magnetic  field  measured 
from  a  bottom-mounted  magnetometer.  Using  a  collocated  pressure  sensor,  and  a  vertical 
magnetic  field  sensor  100  m  away  as  reference  signals,  we  were  able  to  reduce  the  back¬ 
ground  magnetic  noise  on  our  test  signal  by  approximately  18  dB  in  the  frequency  band 
from  5  mHz  to  200  mHz.  The  pressure  reference  signal  was  used  to  reduce  the  coherent 
swell  induced  magnetic  background  noise  in  the  60  mHz  to  200  mHz  band,  and  the  mag¬ 
netic  reference  signal  was  used  to  reduced  the  geomagnetic  background  noise  in  the  5  mHz 
to  60  mHz  band. 

The  magnetic  signature  of  a  transiting  surface  vessel  as  measured  by  a  bottom-mounted 
magnetometer  is  typically  in  the  5  mHz  to  200  mHz  frequency  band.  By  decreasing  the 
noise  on  our  test  channel  by  18  dB  in  this  band,  target  detection  algorithms  may  be  applied 
with  lower  detection  thresholds,  and  target  signals  almost  an  order  of  magnitude  smaller 
than  previously  will  be  able  to  be  detected.  This  will  increase  the  detection  capability  of  the 
surveillance  array  using  this  technique,  and  decrease  false  alarms  produced  by  background 
noise  fluctuations. 

For  surveillance  applications  where  we  have  an  array  of  sensors  or  multi-sensor  platforms 
cabled  together  and  distributed  on  the  ocean  bottom,  the  multi-channel  coherence  technique 
has  great  potential  for  reducing  background  noise  on  the  returned  sensor  data.  For  cabled 
arrays,  all  data  from  all  sensor  may  be  available  simultaneously.  Array  sensors  placed  in 
areas  where  target  vessels  are  unlikely  to  transit  near  could  be  used  as  reference  sensors 
for  the  other  array  sensors  more  likely  to  have  a  target  pass  nearby.  For  surveillance  arrays 
cabled  to  shore,  a  completely  separate  reference  array  could  be  deployed  in  an  area  where 
targets  are  unlikely  to  pass.  This  reference  data  could  be  then  used  to  reduce  noise  on 
sensor  data  collected  from  a  main  array  deployed  in  the  target  transit  area. 

If  a  surveillance  array  has  distributed  multi-sensor  platforms  that  are  not  cabled  together, 
the  multi-channel  coherence  technique  can  still  be  used  to  reduce  noise.  We  have  shown 
that  by  using  a  collocated  pressure  sensor  we  were  able  to  reduce  the  background  noise  on 
a  magnetic  signal.  Other  collocated  sensors,  such  as  three-component  electric  field  sensors, 
could  also  be  used  if  they  show  any  coherence  with  the  background  noise  on  the  test  signal. 

It  should  be  made  clear  that  the  multi-channel  coherence  technique  is  not  just  about  reduc¬ 
ing  background  noise  on  a  magnetic  signal.  We  could  just  as  easily  be  measuring  electric 
fields  or  acoustic  signals,  and  trying  to  reduce  the  background  noise  on  these  signals  by 
using  the  appropriate  reference  sensors. 
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Future  work  needs  to  be  performed  to  determine  the  optimal  sensor  distribution  for  a 
surveillance  array  if  the  multi-channel  noise  removal  technique  were  used.  This  would 
depend  on  the  surveillance  application,  and  in  particular  depend  on  such  things  as;  the 
deployment  depth,  the  target  fields  (magnetic,  electric,  acoustic  etc.)  to  be  sensed,  the  ex¬ 
pected  source  strengths  of  the  targets,  the  different  expected  weather  conditions,  the  tidal 
conditions  and  the  geomagnetic  noise  conditions. 

From  a  computational  perspective,  future  work  could  also  be  performed  in  analyzing  the 
resulting  linear  system  transfer  functions  that  the  model  produces.  Preliminary  study  of 
these  functions  has  shown  that  once  they  have  been  calculated,  they  can  be  used  without 
having  to  update  them  at  every  computation  time  window.  These  functions,  which  in  a 
way  characterize  the  local  environment,  have  some  temporal  stability.  If  this  is  indeed  the 
case,  then  this  will  greatly  reduce  the  amount  of  computation  time  required  to  implement 
the  multi-channel  coherence  technique  in  a  real-time  application.  A  study  of  the  temporal 
stability  of  the  transfer  functions  should  be  performed. 
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